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Abstract: For every Matrix Product State (MPS) one can always construct a 
so-called parent Hamiltonian. This is a local, frustration free, Hamiltonian which 
has the MPS as ground state and is gapped. Whenever that parent Hamiltonian 
has a degenerate ground state (the so-called non-injective case), we construct 
another 'uncle' Hamiltonian which is local and frustration free but gapless, and 
its spectrum is M"^. The construction is obtained by linearly perturbing the 
matrices building up the state in a random direction, and then taking the limit 
where the perturbation goes to zero. For MPS where the parent Hamiltonian 
has a unique ground state (the so-called injective case) we also build such uncle 
Hamiltonian with the same properties in the thermodynamic limit. 



1. Introduction 

It is nowadays widely accepted that the low energy physics of ID locally interact- 
ing quantum systems is captured by the so-called Matrix Product States (MPS) 
[8|[T7] . This idea, rooted in the numerical success of DMRG [30ll23] and based 
on the area-law principle [7] , has had an enormous impact on the understanding 
of ID quantum systems. On the numerical side, it has led to novel methods to 
extract the physics of strongly correlated quantum systems [?7l[251E5LllOlfTT] . On 
the theoretical side, it has served e.g. to understand and quantify the amount 
of entanglement present in ID quantum systems [131 1^15] . or to clarify the role 
of symmetries and its connection with the existence of string order parameters 
[T51ll2j . and ultimately has led to a full classification of ID quantum phases pUl 

A crucial concept in all the above references is the one of parent Hamilto- 
nian, a frustration free gapped local Hamiltonian having the given MPS as its 
ground state. Parent Hamiltonians are extremely useful as they allow us to use 
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the framework of MPS to also study the properties of one-dimensional Haniilto- 
nians. While analyzing systems with discrete symmetry breaking, equivalently 
MPS with a discrete gauge symmetry, it has been observed in [4J, that small 
local perturbations of the MPS do not correspond to small perturbations in its 
parent Hamiltonian, unless the perturbation preserves the gauge symmetry. The 
main aim of the present work is to provide a complete picture of this unexpected 
behavior. Among other things, this could help in the design of MPS-based al- 
gorithms to simulate systems with discrete symmetry breaking in ID or, more 
generally, systems with topological order in 2D. 

In this work, we will study the parent Hamiltonian corresponding to pertur- 
bations of MPS with a gauge symmetry. We will show that in the limit where the 
perturbation goes to zero, the parent Hamiltonian converges to a new Hamilto- 
nian with the following striking properties: (i) It has the same -finite dimensional- 
ground space as the parent Hamiltonian, (ii) it is also frustration free, but (iii) 
it is gapless and its spectrum equals [0, +00) in the thermodynamic limit. These 
properties hold for any MPS with a gauge symmetry, and for generic perturba- 
tions. As opposed to the parent Hamiltonian, we will call this new Hamiltonian 
an uncle Hamiltonian for the MPS. Unlike the parent Hamiltonian, the uncle 
Hamiltonian changes continuously under perturbations of the MPS even in the 
presence of a gauge symmetry. 

As a by-product we obtain new examples which shed new light on the prop- 
erties inherent to critical systems and phase transitions, in the line initiated in 
[ST] . Notice that all our examples are generic, as gapless as they can be since the 
spectrum is the whole positive real line, but they are frustration free and there 
is no power-law decay of correlations within the ground space. Indeed, modi- 
fying our construction, for any injective MPS (which always has exponentially 
decaying correlations), one can obtain simple frustration free models with the 
same spectral properties but with the given MPS as unique ground state in the 
thermodynamic limit - however, the finite dimensional kernels do not coincide 
in this injective case. 

The paper is organized as follows: In Section [2] we start with the basic defini- 
tions and lemmas. Section [3] works out the illustrative example of the GHZ-state 
which helps in understanding the main theorems presented in Section |4j Section 
[5] finally deals with the case of injective MPS. For the sake of clarity, we have 
moved some of the technical proofs to the appendices, where one can also find 
how to treat formally the thermodynamic limit using the GNS construction and 
some basic facts we need about spectral theory of unbounded operators on a 
Hilbert space. 

2. Definitions 

In this section, we will provide the necessary definitions and basic lemmas: We 
start by introducing Matrix Product States, their normal form, and the injectiv- 
ity property. Next, we show how parent Hamiltonians for MPS are constructed. 
We find that parent Hamiltonians change discontinuously under certain pertur- 
bations, which motivates the introduction of uncle Hamiltonians. 



2.1. Matrix Product States. 
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Definition 1 (Matrix Product States). A state |?/;) e (C'')'*-^ is called a 
Matrix Product State (MPS) if it can be written as 



\MiA))= J2 M^n-'-A, 



(1) 



where the Ai, i = 1, . . . , d, are D x D matrices (D is called the bond dimension). 

The matrices {Ai}i can also be thought of as a tensor A with three indices 
iAi)a/3, where the matrix indices a and /3 are referred to as "virtual indices", 
and the index i as "physical index" . Note that this definition of MPS is restricted 
to translationally invariant states with periodic boundary conditions, which we 
will be concerned with in this paper. 

Throughout this paper, we will use a graphical notation for tensor contrac- 
tions such as in Eq. ([T]), cf. [23]: A tensor with k indices will be denoted as a 
box with k legs; e.g., the tensor {Ai)ai3 is denoted as 



a 



1^ 



with the upper leg by convention denoting the physical index. We will generally 
omit the labels. The contraction of two tensors, i.e., summing over a joint index, 
is denoted by concatenating the corresponding indices; e.g., '^p{Ai)ap{Bj)pj is 
written as 



a 



A 



B 



We will also use the more compact notation 

A^B 

instead. In this graphical language, an MPS, Eq. 0, can be expressed as 



A 



Given a tensor T = {Ti)ai3, we will denote by {a\T\l3) the vector 



in the physical space. More generally, matrix operations are generally 
meant to act on the virtual degrees of freedom; e.g., C — A G) B denotes the 
tensor with components Q = Ai (B B^. 

In order to compute expectation values etc., we need to contract tensor net- 
works with their adjoint. Tensors with legs pointing down are always complex 



conjugated!^ 



{Ai)ai:s = {A*)/:sa 



An object of particular interest is the "transfer operator" 



B 



A 



^ Throughout, ■ will denote the complex conjugate and ■* the adjoint operator. 
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which we will occasionally interpret as a map from the right to the left indices, 



2.2. Normal form. MPS are stable under blocking: If we e.g. block pairs of phys- 
ical sites, jk = (i2fe-i,*2fe), the state can again be expressed as an MPS in the 
blocked indices, with tensors Aj^ = Ai^^ -^Ai^^. By blocking a finite number of 
sites P21 and appropriate gauge transformations, any MPS can be brought into 
a standard form [51[T7]. 

Theorem 1 (Standard form for MPS [8',' 17], injectivity). After blocking, 
any MPS can be written in a standard form where the matrices Ai have the 
following properties: 

1. The Ai are block- diagonal: Ai ~ ®f=iA\ (g) Fj, where Aj E (the space of 
Ij X Ij matrices) and the Fj are positive diagonal matrices. 

2. The Ai span the space of block-diagonal matrices: span^Ai = ®f=i -^ij 

3. For all j , the map £j :— E^j has spectral radius one, with 1 as the unique 
eigenvalue of modulus I, and with eigenvectors £j{T} = I and £*(A-'^) ~ A'j^, 
where A\ > 0, tr{A\) = 1. 

Property 2 with every Fj = 1 is called block- injectivity; in particular, if T) = 1, 
and Fi — 1, A is called injective. 

For the rest of the paper, we will always consider MPS in this standard form. 

The uniqueness of the fixed points in each block leads to the following result 
concerning the transfer operators. 

Lemma 1. Let A denote an injective block of an MPS. Then 

k 

I A A \-\ 



A 



+0{e'') 



(2) 



Note that 0(e denotes a bound up to a constant in the exponent. This 
notation will be used throughout this work. 

Lemma 2 (Consequences of injectivity). The following three properties are 
equivalent, 1 <^=> 2 3; 

1. A is injective. 

2. For any X , there exists an \a) such that 



A 



X- 



X 



(3) 



3. There exists a tensor A such that 

^(^^)a/3((A-l),)a'/J' = 
i 

("left inverse to A"). 



A 



- A - 



(4) 
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Also, the following three are equivalent, 4 <^ 5 <^ 6; 
4- ( ^ |j ) is injective 

5. For any X of the dimensions of A, there exists an \a) such that ^ holds, 
and additionally 

® 



- B - 



R 



- L 



, 



and the corresponding statement holdsjor the other three blocks. 
6. There exists a tensor A^^ such that (01) holds, and additionally 



K 



-K 



, 



B 



R 



L 



and the corresponding statement holds for the other three blocks. 

Proof. 1 <^ 2 since by definition, injectivity means that the Ai span the whole 
matrix algebra. 2 3 by defining \aafj) such that in ([3), X — |a)(/3|, and 
choosing {{A~^)i)ap = (aa/jji), and 3 ^ 2 by setting {ai) = tr[{A~^)iX^]. 
4 5 by considering equivalence between 1 and 2 and the matrix 



X 



X 




or the corresponding matrices for the other blocks. 
5 4 since any matrix can be block-decomposed as 

( X 

and for these blocks there exist vectors \ax), \rz)i |W) that give rise to X, 
Y, Z and W when applied to A, B, R and L respectively, and when applied 
to the other blocks. Thus we can consider the sum \ax) + Iby) + \rz} + \lw) to 
satisfy condition 2, and therefore injectivity of tensor in 4. 

5 6 by defining \aap) such that in ([3| X = if both indices correspond 

to the A block or otherwise, and choosing {{A^^)i)a/3 — (aQ^|i), and 6 5 

by setting {a\i) = tr[(A"i)iX^], with X ^ i ^ ^ 



a 



Lemma 3. Under the conditions of Theorem^ the spectral radius p{E'^j) < 1 
fori ^j. 

Proof Let us take X such that Ef^{X) = J^k^i^i^lT = We have then 



that 



\ti:{XA^.X*)\ = \Y,tr{AiX{AirAA,X*)\ = | ^ tr(X(Ai)* v^v^XM^ 



< (^triXiAirA^.AlX*)^ 



1/2 



tr{{AirXA^.X*Ai: 



\tT{XAA^X* 



and therefore for any eigenvalue we have |A| < 1. Note that the inequality is 
strict since we have that span{(/c|yl'|Z)} n span{(m|yl^ = {0} due to block 
injectivity. □ 
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Definition 2 (Span of a tensor, Projector corresponding to a tensor). 

Given a tensor {Ti)^^ with two virtual indices a, /3, and one physical index i 
(which can he a blocked index), we define the span of T as 

span{r} := span { ^ tr[TiX] \i)\X gMd} ■ 

i 

Also, we define the projector corresponding to T, n[T], as the orthogonal pro- 
jector onto spanlT}-"-. 

In particular, T can arise from blocking two or more tensors A from an MPS, 
T = A^A---. E.g., 

span{^-^^} = span { y^^tT[Ai A jX]\i,j)\X e Md} ■ (5) 

Lemma 4 (Gauge transformations for span). Let C : Aio —>■ -Md be an 

invertible map on D x D matrices. Then, span{T} = span{£(T)}, and equally 
n[T] = n[jC{T)], where the natural action of C on three-index tensors {Ti)ap is 
given by [L{T)\i = L{Ti) . 

Lemma 5 (Continuity of projector of a tensor). Let T(e) be a family of 
tensors. If T{0) is injective and T[s) is continuous around 0, then 7I[T(e)] is 
continuous at 0. 

More generally, if T(0) is block-injective and T{e) is continuous and block- 
diagonal in the same basis around 0, then n[T{e)] is continuous around 0. 

Proof. The proof follows directly from the fact that a basis X). of the space of 
(block-)diagonal matrices yields a continuously changing basis of span{r(e)} by 
virtue oiXk ^ J2,tr[T,{e)Xk]\i) = \vk{e)). 

From these bases, the Gram-Schmidt orthogonalization process leads also con- 
tinuously to orthonormal changing bases of span{r(£)}, say {|efc(£))}: Beginning 
from |ui(e)) ^ 0, its normalization is continuous. Then, orthonormalization of 
\v2{s)) with respect to |ei(£)) only involves projecting onto the complement of 
span{|ei(£))} and normalization, and both operations are continuous as long as 
we keep in every moment away from the null vector (this is where the impor- 
tance of injectivity or block-injectivity stems, in the need of the dimension of 
span{T(0)} to bo kept locally). One can also consider the process as projecting 
\v2{e)) onto span{|ei(e))}, and then substracting the result to |w2(e)), which 
yields a non-zero vector that must be normalized. With the rest of the vectors 
the process is similar with finitely many steps, which makes the orthonormal- 
ization process continuous if the injectivity conditions are held and every step 
keeps away enough from 0. Note that injectivity of T(0) and continuity of T at 
is enough to ensure injectivity of T(e) for small s. 

The final projection n[T{e)] is also continuous in e, since projecting any 
vector onto every linear space span{|efc(£))} and substracting the result to the 
vector are continuous in s. 

In the case T(0) is block-injective, the fact that T{e) is also block-diagonal 
makes the rank of {\vk{s))}k constant around 0. This allows to guarantee the 
same result about continuity. □ 
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2.3. The parent Hamiltonian. Let us now turn towards Hamiltonians for MPS. 
We will restrict to translational invariant Hamiltonians, and denote the local 
terms by lowercase letters, e.g., h. When necessary, subscripts indicate the sites h 
acts on, e.g., hi^i+i. We will identify the local operator h and the global operator 
ft, (g) I. The total Hamiltonian will be denoted by the corresponding uppercase 
letter, e.g., H = ^h := Generally, indices will wrap around the 

ends of the chain (e.g., here iV + 1 = 1). 

Every MPS \M{A)) induces Hamiltonians to which it is an exact ground state: 
The reduced state on, say, two adjacent sites, p2, is supported on spanjA— A}, 
and thus for h := n[A^A], h\M{A)) = 0. If D'^ > cP - this can be achieved 
by blocking, and is the case for the standard form of Theorem fll - n[A—A\ 
is non-trivial, and we obtain a non-trivial Hamiltonian with two-body terms h 
which has \M{A)) as its ground state. 

If we use \M{A)) in the standard form of Thm.[l| this Hamiltonian is partic- 
ularly well-behaved: 

Definition 3 (Parent Hamiltonian). Let \M{A)) he a (block- Jinjective MPS, 

11 and let h = n[A—A]. Then, the Hamil- 
lamiltonian. 



i.e., satisfying condition 2 of Theorem 
tonian H = is called the parent . 



Theorem 2 (Ground state space of parent Hamiltonian |15LI17L[24] ). 

The parent Hamiltonian has a T>-fold degenerate ground state space spanned by 
\M(A^)); in particular, \M{A)) is one of its ground states. Also, the parent 
Hamiltonian is gapped in the thermodynamic limit. 

Remark 1. Note that in order to ensure the correct ground state subspace in The- 
orem[2] a weaker condition than the injectivity of each tensor used in Definition|3] 
is enough: It is sufficient to take projectors onto the orthogonal complement of 
the span of fc -I- 1 sites, where k is chosen such that blocking k sites makes the 
tensor injective |22) . 



2.4. Uncle Hamiltonians. The parent Hamiltonian construction can be inter- 
preted as a map from the set of MPS to the set of Hamiltonians, H : A^ H{A), 
which associates to any MPS \AI{A)) its parent Hamiltonian H{A). While this 
map is well-behaved in terms of the properties of H{A), we are also interested 
in its continuity: If we change A smoothly, A ^ A + eP, does H{A) change 
smoothly as well? If it were so, this would allow us to study perturbations of 
the system by looking at perturbations of the MPS tensor A. For injective MPS 
in their standard form, Lemma [5] tells us that this is indeed the case. On the 
other hand, if A is block-injective. Lemma [5] requires P to be block-diagonal as 
well, and it is indeed easy to see that random perturbations P will lead to a 
discontinuous change in H{A + eP), as rank(span{^ — A}) is different for e 7^ 0. 

This discontinuity motivates the introduction of uncle Hamiltonians, which 
are robust under specific perturbation of the MPS tensor. 

Definition 4 (Uncle Hamiltonians). Let \M{A)) be a MPS in standard form, 
and let P Md- Then, the imcle Hamiltonian induced by P is the Hamiltonian 



H'p := lim H(A + eP) , 

£-i.O 



(6) 
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whenever this limit exists. This is, the uncle Hamiltonian is the limit of the 
parent Hamiltonian of the perturbed MPS for the perturbation going to zero. 

For A an injective tensor, the limit exists and is equal to the parent Hamiltonian, 
following Lemma [5j Parent and uncle Hamiltonian coincide for injective MPS. 
Thus, we focus our attention on the uncle Hamiltonians of block-injective MPS; 
most of this paper is devoted to studying their properties. 

In the following, H — ^h will denote the parent Hamiltonian, while H'p = 
^ h'p denotes the uncle Hamiltonian, where we will occasionally omit the sub- 
script P. 

3. Example: The GHZ state, the Ising model, and the XY model 

We will start our discussion of the properties of the uncle with the GHZ state, 
which has the Ising model as its parent Hamiltonian. The unnormalized GHZ 
state can be expressed as an MPS 

|GHZ) = |00---0) + |ll---l),= tr[A,, ...A,;J|ii,...,iL) 

n,...,iL 

with ij £ {0, 1}, where Aq = ( J []) and Ai = Following Definition [3| the 

parent Hamiltonian of the Ising model can be constructed from the span of two 
sites, 

span{A^A} = span{|00), |11)} , 

which is indeed (up to an additive constant) the well-known Ising Hamiltonian 
il-[|00)(00| + |ll)(ll|]. 

3.1. Uncle Hamiltonian for the GHZ state. Let us now construct the uncle 
Hamiltonian for the GHZ state. According to the definition, we first need to 
fix a perturbation P of the MPS tensor A, 

Next, we need to consider the MPS \M{Ce^p)) with C^^p = A+eP, and construct 
its parent Hamiltonian. For a generic P, we need to block two sites to reach 
injectivity. Thus, following Remark [T] we need to construct the terms of the 
parent Hamiltonian as the projector onto the complement of the span on three 
sites, 

53(e) = span{Ce,P^Ce,P^C,,p} . 

S^{e) is spanned by the four vectors Vap = (a|C'e,p — Cg^p — C^^pI/?), which are 
straightforwardly found to be 

VQo = |000) + 0(e) , 

vm = e[&o|000) + (6o + &i)|001) + % + 6i)|011) + + ©(e^) , 

vio = e[co|000) + (co + ci)|100) + (cq + ci)|110) + ci|lll)] + 0{e^) ; 
vii = |111) +0(e) . 
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If 60 + ^1 7^ and cq + ci 7^ (which happens almost sureljQ), this can be 
transformed into an ahernative set spanning S^^e), 

{|000) + 0(e), |001) + |011) + 0{e), |100) + |110) + 0{e), |111) + 0{e)} . (8) 

The parent Hamiltonian of the perturbed MPS |M(Ce^p)) is thus i/p,^ = J2 hp,e, 
with each hp^ acting locally on three consecutive sites, and projecting onto 

In order to obtain the uncle Hamiltonian we finally need to take the limit 
e — >■ 0. Then, the four states in Eq. Q become orthogonal, and the family hpj, 
converges to the projection onto the orthogonal complement of 

span{|000), |0+1), |l+0), |111>}. 

Here, |0+1> = |0>|+)|1>, with |+) = (|0> + |l))/\/2. Thus, the uncle Hamiltonian 
has local terms 

h'p = l- [|000)(000| + |111)(111| + |0+1)(0+1| + |1+0)(1+0|] (9) 

Note that this limit does not depend on the perturbation P (as long as feo+^i 7^ 
and Co + ci 7^ 0), and will be called h' or h[_i ^ ^^i in the following. The parent 
Hamiltonian W' is obtained as the sum H' =^h' = ^ ^j^^. 



3.2. Ground space of the uncle Hamiltonian. What is the ground space of the 
uncle Hamiltonian? Since 

ker(/i') = span{|000), |111>, |0 + 1), |1 + 0)} D span{|000), |111)} = ker/i , 

the ground states |0 . . . 0) and |1 . . . 1) of the GHZ parent Hamiltonian are also 
ground states of the uncle; in particular, the uncle is frustration free. However, 
h' allows for additional ground states. Indeed, if we consider the ground state 
space of h' acting on m consecutive sites (with open boundaries), the ground 
space is 

span{|0...0),|l...l),^|0...01...1),^|1...10...0)} ^(€2^", (10) 

where the sums run over all positions of the "boundary wall" 01 and 10, respec- 
tively. Yet, when closing the boundaries, the additional states ^ |0 • • • 01 • • • 1) 
and ^ |1 • • • 10 • • • 0) stop being in the intersection of the kernels, and the ground 
space of the uncle Hamiltonian coincides with the ground space of the parent 
Hamiltonian. Intuitively, with periodic boundary conditions the boundary walls 
need to come in pairs, and it is impossible to give both of them momentum zero 
as they meet which is not in the ground space of h' . 

^ Note that whenever we say something happens almost surely means not only that it 
happens with probability 1 but also that the set of perturbations which may not satisfy the 
statement form a closed algebraic variety of dimension strictly lower than the set of all possible 
perturbations. 
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3. 3. Spectrum of the uncle Hamiltonian. Let us now show that the uncle Haniil- 
tonian of the GHZ state is gapless. To this end, we consider the unnormaHzed 
states 

-A'<i<-1 
l<i<Af 

with 

= |o-^-io---o'i---io^ •••0^+1) 

on a chain of length 2N + 3, where the superscripts indicate the position of the 
corresponding qubit. (This notation will be used throughout this work.) These 
states are orthogonal to the ground space, and (^atI^at) = iV^. Further, \(t>N) is 
almost a zero momentum state of the two boundaries: It is a ground state of all 
terms in H' except /i_i,o,i, /iAf,7v+i,-A'-i, and /iat+i.-tv-i.-jv, and by counting 
the violating configurations, we find that ((/)Ar|iJ'|(/)jv) = 0{N). Hence, for the 
energies of these states we have 

iMH'\M ^0{1/N). (12) 



This implies that (on a chain of length 2N + 3) H' has at least one eigenvalue 
'^N < ^'^(yin JwT^ ^ '^(1/^)' i-^-i family of uncle Hamiltonians is gapless. 

Does H' have a continuous spectrum? One idea to prove so would be to 
to construct momentum eigenstates with an energy scaling like 0{k'^ /N"^) {0 
denotes the exact scaling rather than an upper bound) . To this end, we give the 



sum in Eq. (11 1 a momentum, 

J2 e^.^fc/A^i^^^^^ (^3) 



ON,k 



-N<m<-1 
l<n<N 



It is straightforward to see that in addition to the 0(N) contribution from be- 
fore, the — 1 terms h—N — 1, —TV, —N -I- 1, ... , h-2,-1,0 all give a contribution 

0{k^N), resulting in an energy ^^JnIvn)^ = ^(k'^/^^)- 

Unfortunately, the existence of states with energy 0{k'^ /N"^) does not allow 
to conclude that the spectrum of H' is dense: The existence of a state with 
energy E only implies the existence of eigenvalues Xi < E and A2 > E, but 
tells us nothing about their exact value. (The reason this worked for the gapless 
excitations was that H' > 0, and that |0jv) was orthogonal to the ground space.) 
It is, however, indeed possible to show that H' has a continuous spectrum. For 
the GHZ example, this can be done by mapping it to the XY model as discussed 



in Sec. 3.5 and in Sec. 4.4 we give a proof that uncle Hamiltonians of arbitrary 
MPS have a continuous spectrum which works directly in the thermodynamic 
limit. 



3.4- Gapless uncles for unique ground states. Can we obtain uncle Hamiltoni- 
ans with similar properties in the case of MPS which are unique ground states? 
Lemma [5] tells us that this cannot happen as long as the MPS tensors are injec- 
tive, which is always the case as long as such an MPS is in its standard form: In 
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that case, the uncle Hamihonian is equal to the parent Hamiltonian. However, 
as we will demonstrate in the following, interesting uncle Hamiltonians can be 
obtained by choosing a different MPS representation. 

Consider a qubit chain (C^)**^, and a state \M{A)) = |0. . .0). Clearly, this 
is a unique ground state of a gapped local Hamiltonian, with standard MPS 
representation Aq = (1), Ax — (0). However, we can write the same state with 
bond dimension 2 and 

^"=(01) ^'^^^i=(oo) 

We can now perturb A with a perturbation P as in Eq. ([t]), with Oi = di and 
bi = Ci, i = 0,1. A calculation similar to the one for the GHZ state shows 
that in the limit of a vanishing perturbation, the ground space on three sites 
is 53 = span{|000), \W)), where \W) = (|001) + |010) + |001))/%/3 (as long as 
61 7^ or ci 7^ 0); the uncle Hamiltonian h' is the projector onto , 

h' = 1- [|000)(000| + . (14) 

For an open chain of length H' ~^h' has the two ground states 

|0„) = I00---0) and 
|W„) = |100 • • • 0) + |010 • • • 0) + . . . + |000 • • • 01) . 

Different from the GHZ case, the extra state \Wn) does not disappear from 
the kernel when closing the boundaries - the uncle Hamiltonian on a chain 
with periodic boundaries has a two-dimensional ground space span{|OAr), iWOv)}. 
Note, however, that the thermodynamic limit of \Qn) and \Wisi) is the same, and 
thus, the ground space collapses to the original one in the thermodynamic limit. 
Again, we can construct gapless excitations by considering the states 

|0„) = ^ |0 • • • orOO • • • OPO • • • 0) . (15) 

As before, they are orthogonal to the ground space, and their energy is 0{1/N). 
Alternatively, we could have also choosen a W state with momentum, 

Iv'Jv.fe) -^e^^lO-.-OOl^OO-.-O), fc^O , 

j 



which are also orthogonal to the ground space, and have energy 0{k'^ /N"^). 

Again, the spectrum is dense in the thermodynamic limit. This can once more 
be verified by a mapping to the XY model or, directly in the thermodynamic 
limit, using the methods described in Sec. 4.4 
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3.5. Relation to the XY model. Both the uncle Hamihonian for the GHZ state 
and for the |0 . . . 0) state are closely related to the XY model (or, equivalently, to 
non- interacting fermions) , which can be used to immediately infer that they are 
gapless models with continuous spectra. Let us first consider the uncle Haniilto- 
nian of the GHZ state, Eq. (|9]): It can be rewritten as 

h' = -\[l(E) Z (g) Z + Z (g) Z (E)1 + I(E) X (E)l - Z (g) X (g) Z] + ll(E)I(E)I . (16) 

This is exactly the Hamiltonian discussed in Eq. (11) from ^31j at <? = (|^ It can 
be solved either by transforming it to non-interacting fermions, or by a duality 
transformation to the XY model [T^ |^ The resulting Hamiltonian is 

HxY ^-l^[X^g) X^+i +Y,g) Yi+i + 2Z,]+ const. 



Indeed, this point in the XY model, which can be solved exactly by mapping 
it to non-interacting fermions |14j . is known to be gapless with a continuous 
spectrum. 



Let us now turn to the uncle Hamiltonian ( 14) for the |0 . . . 0) state. Let us 
first replace the uncle Hamiltonian with a simpler one with the same spectral 
properties. Namely, let 

/i' = 1- |00)(00| - |^+)(<Z>+| , (17) 
with \<P+) = (|01) + \10))/V2. We have that 

which implies that for any finite chain, the ordered eigenvalues of iJ' = ^ h' 
and Xi oi H' — J2 ^' ^''^ related by 

■^Xi < \i < Xi . 

I.e., if we want to determine essential spectral properties of H' , such as whether 
it has a continuous spectrum, we can equally well study H' . Since h' can be 
rewritten as 

h' = -l[Xg)X + Yg)Y + Zg)l + lg)Z] + ^Ig)l, 



this yet again gives rise to the same point of the XY model, Eq. (16), proving 



that the Hamiltonians Eqs. (14) and (171 have a continuous spectrum 



^ In fact, the construction in Eq. (10) of j31| is, up to a gauge transformation, equivalent to 
the uncle construction, with e = ^/g. 

■* The partial isometry T : |2i , . . . , ijv) l*i + *2, ■ ■ . , *]V + h) from the eigenspace oi 
associated to the value 1 to the even Z parity space induces the duality mapping Xi i— > 
ToXiO T-l = (gi Xi and Zi (gi Zj+i i-s> T o (Z^ (gi Zi+i) o T'^ = Zi. 
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4. Properties of the uncle Hamiltonian 

In this section, we will see that the observations we made for the GHZ uncle 
Hamiltonian generalize to uncle Hamiltonians of arbitrary MPS with degenerate 
ground states (under some generic conditions): Their ground state space is equal 
to the ground state space of the parent Hamiltonian, they are gapless, and they 
have a continuous spectrum. This section will also contain the proofs which have 
been omitted for the special case of the GHZ state. 

For simplicity, we will focus here on the case where the MPS tensor Ci in its 
standard form, Theorem[l] has two blocks, Ci = Ai(BBi, but the same procedure 
can be followed in the general case: The results are completely analogous in case 
of multiple different blocks, but there are some differerences if there are blocks 
with a multiplicity larger than one. We will comment on this particular case in 
Section [Sl 

Thus, in this section we will be dealing with an MPS |A/(C)), 

where both A and B are injective. We will choose A and B in the normal form 
of Theorem [T] The parent Hamiltonian of this MPS consists of local projectors 
n[C^C] with kernels 

span{C^C} = spa.n{A^A} + span{B^S} , 

where the two-dimensional ground state space is spanned by \M (A)) and \M{B)) . 



4-.1. Form of the uncle Hamiltonian. We will now determine the general form of 
the uncle Hamiltonian (including whether the limit in its definition exists). 

Theorem 3. Let \M{C)), be an MPS with block- injective C — A Q) B , and let 

be an arbitrary tensor, such that the "uncle tensor" 

( A^A A^R+R^B\ 
^ " \B^L + L^A B^B ) ■ ^^^> 

is injective. Then, the uncle Hamiltonian induced by P exists and is a sum of 
local terms hp — n[U]. 

Proof Consider a perturbation = C + eP of the MPS |M(C)). We have that 

rjLr^f A^A + Oie) e{A^R + R^B) + 0{e^)\ 

~ \e{B^L + L^A) + 0{e^) B^B + 0{e) J' ^^^> 

clearly, for e 7^ the map which multiplies the off-diagonal blocks by 1/e is 
invertible and thus (following Lemma |4]) 

n[Ce^Ce] = n[c,{c,^Ce)] - n[u + o(e)] . 

Following Lemma[5j the limit lime_>.o n[U + 0{e)] exists whenever U is injective, 
and equals n[U]. □ 
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The required injectivity of the tensor in Eq. (181 follows in particular from 
the following condition on the perturbation. 

Definition 5 (Injective perturbation). A perturbation P of an MPS \M{C)) 
(in the notation of TheoremW^ is called injective if 



A R 
L B 



(20) 



is an injective tensor. 

Lemma 6. // a perturbation P is injective, then the resulting "uncle tensor" U , 



Eq. (18), is injective. 



Proof. Let us consider condition 5 from Lemma [2] and any X of the dimensions 
of A— A. Since the perturbation tensor is injective, there exist vectors \a) and 
\a') such that 

® ® ® ® 



A 



= I , and 



B 



R 



L 



= 



The product \a) ® \a') yields then for the uncle tensor 



®© 



- u 



X 




A similar reasoning can be followed for the rest of blocks, thus satisfying 
condition 5 from Lemma [21 □ 



A perturbation P is generically injective if d > D^, and we will assume 
injective perturbations in the following. Note that unlike for the |GHZ) state, 
the uncle Hamiltonian does in general depend on the perturbation (though only 
on its off-diagonal blocks R and L)^ 



^ If D < d < and U [Eq. | |l8| ] is injective one can construct both parent and uncle 
Hamiltonians from projectors onto the span of three consecutive sites - this is what we did for 
the GHZ example. Then, the uncle Hamiltonian is the projector associated to 



- B^L^ 



L^A^A 



A^R^ 
B^B^ 



B - 
B 



RSLB^B~- 



While we will restrict to injective perturbations for clarity, the same steps can be followed 
assuming only injectivity of U . 
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4-2. Ground space of the uncle Hamiltonian. In the following, we study the 
ground state space of the uncle Hamiltonian. Throughout, we will restrict to 
injective perturbations. 

Theorem 4. Let P be an injective perturbation of an MPS \M{C)), C — AOiB. 
Then, the ground space of the uncle Hamiltonian Hp is spanned by \AI(A)) and 
\M{B)), and thus equal to the ground space of the parent Hamiltonian. 

Proof. The parent Hamiltonian is frustration free, i.e., its ground states minimize 
the energy of each local term. The ground space is thus 



ker(iJ) = P|ker(/i) . 



Since ker(/i) C ker(/ip), it follows that 



kei{H'p) = P|ker(/i^) D ker(ff) , 

i.e., the uncle Hamiltonian is frustration free, and any ground state of the parent 
is also a ground state of the uncle. 

In order to classify all states in ker(_ffp) = p|ker(/ip), we will follow the 
same steps as for the proof of the ground space structure of the parent Hamilto- 
nian j24] : First, we will prove inductively how the ground space on a chain with 
open boundaries, P| *L ker /ip^ grows - the intersection property. Then, we 
will show how the ground space changes when we close the boundaries - the 
closure property. 

Lemma 7 (Intersection property). Given a chain of length n, let Sn be the 

vector space 



Sn 
Bn 



Ar 



Bn + Rn 



Ln , where 



X 



^ /XeMi}, 



B - - B 



X 



Ln — 



pos R 
pos L 





- R - 


- B 


X 


B 









B. 


L 


- A 


X 


A 









IX e A^;xm}, 

IX e Almx/}, 



(21) 



where the sums run over all possible positions of the R or L, padded with A 's 
and B ' to the left and right as indicated. Further, let ( ^ p ) be injective. Then, 
the intersection property S'„ (g) n (8) S"™ ~ Sn+i holds. 
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The proof of the Lemma is given in Appendix [A) 
Starting from ^2 — ker(/i']^ 2): ^-nd using that 

k k-1 k 

f|ker(/.^^+i) = fl ker(/j:_,+i) n f]kev{h'^^^^,) , 

1=1 i=l 1=2 

the lemma allows to inductively prove that the ground space on n consecutive 
sites with open boundaries is HlLi^ ker(ft,^ j^j^) = 5„: It differs from the ground 
space of the parent Hamiltonian by the presence of the "zero momentum domain 
wall states " i?„ and L„ , analogous to the domain wall states for the GHZ uncle. 
It remains to show that these states disappear from the kernel when closing the 
boundaries. 

Lemma 8 (Closure property). Consider a chain of length N, and let S\cit = 
Sn defined on sites 1, . . . , N , and S'right — Sn defined on sites 2, . . . , iV, 1, using 
the definitions of the previous lemma. (I.e., for S'right the ordering of sites is 
shifted cyclically by one.) Then, 

^icft n bright = span{|A/(A)), \M{B))}. 

The proof is again given in the Appendix [Xj 

The closure property shows that if we close the boundaries on a chain of 
length N, we indeed recover the ground space of the parent Hamiltonian, since 

N N-1 N 

n kcr(/^- ,+1) - fl ker(/i^_,+i) n f| ker(/i^^,+i) = Sieft n Snght • 

i=l i=l 1=2 

Together, the two lemmas thus prove Theorem |4j □ 



4.3. Gaplessness of the uncle Hamiltonian. One of the key properties of the 
parent Hamiltonian is that it exhibits a spectral gap above the ground space 
jl5j . On the other hand, as we will prove in the following the uncle Hamiltonian 
is generically gap less: 

Theorem 5. The uncle Hamiltonian Hp is gapless for almost every P. 
The Theorem can be proven using the following Lemma: 
Lemma 9. For a chain of length 6N +1, let 



-2N<i<-N 
N<j<2N 



where 



(22) 



\C^, 



r A 



- R - B 



B - 



- A 



A 



(23) 



Then, for almost every R and L ( and thus almost every P ), the following holds: 
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1- {^nIM =OiN^). 

2. {M{A)\<j>N) = 0(e-^) and {M{B)\(Pn) = 0{e-^). 

3. {MH'p\(i3N) ^0{N). 

Here, 0{-) denotes both lower and upper bounds on the scaling. 

Note that \4>n) generalizes the GHZ "boundary wall" ansatz for low energy 



states, Eq. ( 11 ). The range for i and j in ( 22 ) is chosen such that R and L move 
over a region of size N each, leaving two separating regions of length 27V each 
which contain only A or B tensors, respectively. 



Proof (of Theorem^. For the normalized states \4)n) '■— \4>n) /\\\4'n)\\ on a 
chain of length QN + 1 (with |(/)jv) of Lemma[9|, we have that \(f)M) tends to be 
orthogonal to the ground space of H'p and {4>n\H'p\4)n) — > as — ?> oo. To- 
gether with simple spectral decomposition arguments, this implies the existence 
of a sequence Jat — !• such that H'p (on 6A^ + 1 sites) has at least one eigenvalue 
in the interval (0, 5n)- □ 

Proof ( of LemmaW^. 

1) {4>n\4>n) — 0{N^): Results from Lemma [l] allow us to approximate 



(CijICfe.i) 



A- - A- R-, r 
B I 1 R U R ur I- 



B - B 



B 



B 



B 



L - A 



A- A 



0(e 



-2N\ 



(24) 

for i < k and j < I, and correspondingly for the other cases. [Note that we have 
used p(Eq) — p{E^) < 1, Lemmajsj together with Eq. ^ to bound the error 
term.] It follows that {4>n\4'n) = ^r^l + 0(e~^^), 



with 



- A 




R 




1 

- B 




1 

B 





T-| A[- \n=0 



- B 



B 



- A 



R ^ 



and correspondingly for Using 



N 



NI- 



B 



we find that = CrN + 0(1), with 



m 



0{l) 



R 






A 




R 






R 






B 












1 










1 


R 






R 




B 






A 






R 



Cp is a quadratic function in R which does not vanish identically (e.g., there 
exists an R for which HrHa = and the first term is non-zero). Thus, the R 
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for which Cr = form an algebraic variety of smaller dimension, and Cr 7^ 
for almost all 

The same argmnent can be used to see that Sl = ClN + 0{1) where Cl 7^ 
for almost all L, and thus, {4>N\(f>N) — 0{N'^) for almost every perturbation as 
claimed. 



2) {M{A)\c^n) 
{M(A)\M 



0{e-^) and {M{B)\(I)n) = 0(e~^); In the scalar product 



E 



A 



- A 



A 



A 



- A- A 



1 



-2N<i<-N 
N<j<2N 



R 



B 



B 



A 



(25) 



every summand contains (Eg)"^-^ . Using Lemma and the fact that there are 
only 0{N^) summands, {M{A)\(I)n) — 0(e^^) follows, and analogously for 

3) {(t)]si\H'p\4)iq) ~ 0{N): The only terms in H'p which give a non-zero energy 



are h' 



2N-1, 



h' h' 

-2N' '^-N,-N+l' "■N-1,N 



, and hy 



2N.2N+1' 



For each of them, N 



summands in (22) contribute, and thus, {4>n\H'p\4>n) = 0{N). 



- 1 
□ 



The spectrum of the uncle Hamiltonians is M^. In order to study more 
properties of the spectra of the uncle Hamiltonians we need to move on to 
the thermodynamic limit. A formal description of it via GNS-representations 
with respect to ground states can be found in Appendix [Bj The spectrum in 
the thermodynamic limit can be found to be the whole positive real line and 
the spectra of the finite sized chains can be proven to tend to be dense in the 
positive real line. 

Through the GNS-representation the problem is translated into studying the 
action of H' on the space 



S 



i<j 



where S, 



A- 


- A- 


X 


-A 


- A 



X 



(26) 



j+i 



and the spectrum of its unique self-adjoint extension 

UJA 



S 



S, with uj = 



A- 



A 



the ground state used for the representation. There 

exists a unique self-adjoint extension because S* is a dense set of analytic vectors 
for H' , and therefore H'\s is essentially self-adjoint [51] . 

In first place, we must show that H'^ is gapless. A family of states related 
to those we used previously for finite chains in Lemma [9] and Theorem [s] let us 
show the absence of gap: 



E 

-2N<i<-N 
N<j<2N 



\C^, 



(27) 



^ If TTa^s = 0, i.e., Eg = 0, such as for the GHZ state, one can prove that Cu ^ for any 
injective perturbation. 
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where 



A- A- R- B 



B - L - A- A 



(28) 



Theorem 6. The uncle Hamiltonian Hp is gapless in the thermodynamic limit 
for almost every P. 

Proof. We consider the operator H'^, which models the thermodynamic hmit of 
Hp with boundary conditions described by the A tensor. 

The bounds from Lemma [9] also apply to the states in (27 1, and use of the 
spectral theorem for unbounded operators (Appendix [C| lead to the fact that 
H'^ is gapless. □ 



This last result shows the existence of a sequence of elements in a{H'^) tending 
to . They will allow us to prove, together with the next two results, that the 
spectrum is the entire positive real line. 

Proposition 1. A real value A € spec(i?^) iff there exists a sequence of nor- 
malized states {\(p\,k)}k G S such that \\{H'^ — A/)(|(/3A,fe))|| ~^ 0. 

Proof. This follows from the fact that the residual spectrum of i/^ is empty 
-since it is self-adjoint-, the definition of a value lying in the point or continuous 
spectrum (Appendix [C |, and the f act that H'^ is the closure of H' acting on S, 
that is, graph(ij;;,) graph(iJ'|5) C S x S. □ 

Proposition 2. // some values {Xi, . . . , A„} lie in the spectrum of H^ then the 
sum Xi also lies in the spectrum of H'^ . 

Proof. We will state the proof for two values a = Ai and 6 = A2. The construction 
for the sum of more values is completely analogous. 

For both a and b we can find some sequences of normal states {Iva.fcljfe and 
{\^b,k)}k C S verifying the previous proposition, with \\H[^{\ipc,k)) - \^c,k)\\ < 
1/k for c — a,b. We can assume, due to translationally invariance of H' , that the 
first sequence lies in Ui<-fc ^i,-k, and the second one is contained in ljj>fc ^k,j- 
These states would then have the form 



Wa.k) 



A- 


- A- 




- A- 


- A 



m,k) 



A- 


- A- 




- A- 


- A 



for some tensors Xa,k, ^b.fc- The normalization conditions would be 



(M I i- --'T-r-^ (T) = 1 = (a) I I - -- T-r-^ rp 

Xa,k Xb,k 
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From these states we can consider the 'concatenated' states 



A- 


- A- 




- A- 


- A 




A- 


- A- 


Xb,k 


- A- 


- A 



(29) 

Note that the separation between the X blocks is increasingly growing. 

And due to the structure of H[j the image of lJi<-fe ^i-k is contained in 
Ui<-fe 5'i,-fc+i, and i?^(Uj>fe Skj) C \J.^^ Sk-i,j. Moreover, there exist tensors 
X'^ f. and X'^ f, such that 



H' 



A- A 



I I I I 



A- A 



aH 



■k+2 



A- 


- A- 


Xb,k 


- A- 


- A 


- u 


A- 


X'b,k 


- A 



These new tensors also allow us to describe the image of the concatenations: 



A- A 



I I I I 



X, 



a,k 



A- A 



A- 



_I_L 



_I_L 



X' 



a,k 



- A 



-A 


A 


Xb,k 


- A 


- A 


k-1 

1 1 1 1 1 1 1 1 


A- 


A- 


Xb,k 


A- 


A 



I I I I 



A. A. -Xa k 



k-l 

-U I I 



A- 



I I ^ 



X', 



b,k 



- A-- 



k-2 



Let us call 1^'^. „) and I'P'i^ h) these two summands. 

We then have that \\HUm) - (a + b)\^k)\\ < WW^^J " a|<^fe)|| + 

We can derive a bound for the first summand: 



A- 


X'a,k 

1 1 1 1 1 1 


- A 


1 

A- 


X'a,k 


1 

- A 



A- A 



Xb,k 
1 1 1 1 


- A- 


- A 


Xb,k 


1 

- A- 


1 

- A 



+|a|2(l+0(e-*=))- 



-2 Re 



A- 


A- 


Xa,k 


- A- 


- A 


1 

A- 


1 II II 1 

X'a,k 


1 

- A 



A- 

1 


- A- 

1 


Xb,k 
1 1 1 1 


- A- 
1 


- A 

1 


A- 


- A- 


Xb,k 


- A- 


- A 



Xa,k Xbjt 
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-2 Re 



A- 


1 1 1 1 


- A 


X'a,k 



X, 



TIE 



b,k 



ITX 



+ 0(e-'=) = 



Xb,k 



X' 



a,k 



1ZE 



|ar-2Re 



X'a.k ^ 



A- 


Xa^k 
1 1 1 1 


- A 


X'a,k 




0{e-^) < l/e + Oie-"), 



where Re(-) denotes the real part. 

A similar bound can be found for the second summand, and we can derive 
the bound 

\\K{\^,))-{a + b)\<P,)\\^0{l/k). 

We also have that ^ 1- Therefore, the sequence |^fe)/|| || satisfy 

the conditions in Proposition [l] for a + = Ai + A2, and consequently this sum 
lies in the spectrum of H^. 

Longer concatenations would prove the result for any finite sum among values 
from {Ai, . . . , A„}. Note that the bound we get depends on the number of values 
Xi we are summing, 



\Ki\<P,))-iJ2^^)m\\=Oin/k). 



□ 

Theorem 7. The spectrum of the uncle Hamiltonian Hp in the thermodynamic 
limit is the whole positive real line for almost every perturbation P. 

Proof. The set of finite sums of any sequence of real numbers tending to is 
dense in M+. Since there exists a sequence of elements in (j{H[J) tending to -as 
it can be deduced from Theorem[6]- and any finite sum of these elements lies also 
in a{H'^), which is closed, this last spectrum must be equal to M"*". Therefore 
the spectrum of Hp is the whole positive real line. □ 



Spectra for finite chains. After this discussion on thermodynamic limit Hamil- 
tonians we need to go back to the finite chains, and study how the spectra of 
the uncle Hamiltonians constructed on finite spin chains tend to be dense in IR+ 
as the size of the chain grows. 

Given i < j, Sij can be easily mapped to any finite chain of lenght 2A^ -f 1 
for N > max{j,j} via 



Cat : 



X 



j+1 



A 



I I I I 



X 



A 



and this family of maps capture important information since they tend to be 
isometric embeddings. 



A 



Lemma 10. For fixed i and j , cn is an isometry up to a correction 0{e ) 
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Proof. It follows from Lemma [T] 



□ 



Note that in the case of G-isometric MPS 24J, these maps are isometries. 

Since we will need to keep track of how close some vector evidencing the 
existence of elements in the spectrum is from being an eigenvector we introduce 
the following definition. 

Definition 6. A normalized vector will be called an approximated eigenvector for 
an operator A and a given value X, and for an error e if \\{A — XI)\(px^k)\\ < £• 
In the case of non-normalized vectors they must satisfy 



\{A^XI){\^,^k))\\ 

III^A,fe>ll 



< e 



Lemma 11. For any given values X e , n G N and S > 0, there exists a 
value Nq such that we can find approximated eigenvectors for the values jX, 
j = 1, . . . ,n, and for an error at most S for every finite chain with length greater 
that 2iVo + 1. 

Proof. Let us take S' = 5 /An. For A and an error 6' , a normalized approximated 

can be found in some 



A- 


A- 


X 


- A- 


- A 



eigenvector \(p\,s) = 

-M-l M+1 

S~M,M such that \\{H'^ — AI)|(/3a,5)|| < S' . We can now find a value r such that 
the following vectors are respectively approximated eigenvectors, not necessarily 
normalized, for the values jX, j = 2, . . . , n, and an error at most jS' < S/2: 



_1_L 



X 



A 



J_L 



_L± 



A- X 



_u_ 



X 



A 



A 



J_L 



_L± 



X 



A 



_1_L 



X 



A 



_1_L 



X 



A 



_LL 



X 



N-3 



X 



A 



where the r denotes how many A tensors are missing in the diagram, and the 
— 3 refers to the number of X blocks with r A tensors between every two of 

them are also missing. 

A value M' can be found such that all these states belong to S-m',m'- And, 

due to Lemma flO) there exists a value A^'o such that the maps 



Cat 



5_ 



M'.M' 



A- 


- A- 


Y 


- A- 


- A 



5*2^+1 

I I ['■■■■■I I I 
A- Y - A 



make each e{\ip^^\)) an approximated eigenvector for jX and an error ^, for the 
uncle Hamiltonian for the corresponding finite size chain, and for every A'^ > Nq. 
The corresponding normalized vectors are approximated eigenvectors for the 
same values, and therefore satisfy the statement in the lemma. □ 



A 



However, these jA need not be in the spectrum, but indicate the existence of 
elements in the spectrum close to them, as it is shown in the following lemma. 
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Lemma 12. Let A be a self-adjoint operator on a finite dimensional Hilhert 
space, and A a positive real value such that there exists a unitary vector \^p\ s) 
with\\{A-\l)\ipxM<5- T/ien ct(A) n (A-(5,A + (5) 7^ 0. 

Proof. Let {Ai} be the set of eigenvalues (possibly repeated) of A, and a 
corresponding orthonormal basis of eigenvectors. Then \(p\^s) can be written as: 

|<PA,5)=5]a.|<^^), (30) 

i 

and A{\lpx,&)) = 

If cr{A) n (A — (5, A + (5) were the empty set, we would have that 

\\{A-Xl)\ipxM = \\Y,a,{K-\)\m > ||^a,<S|0,)|| =5, (31) 

i i 

which contradicts the conditions of the lemma. □ 

Theorem 8. The spectra of the uncle Hamiltonians for finite size chains tend 
to be dense in the positive real line. 

Proof. For any given values L,m e N, we can set n = Lm, A = 1/m and 
6 — l/{3mN) — l/{3Lm^). For these values, approximated eigenvectors can be 



found following Lemma 11 which indicate, due to Lemma 12 that (jX — S,jX - 
S) n <j{H') ^ 0, j = 1, . . . , n, for every long enough chain. Therefore eigenvalues 
for every long enough chain can be found distributed in disjoint intervals centered 
on the set j/m, = 1, . . . , mL as depicted in the following diagram. These sets of 
eigenvalues, however possibly different for every chain length, tend to be dense 
in M"*", as we consider higher values for L and m. □ 



1/m 

I I I I I I I I I I I I I I I I I 



2 3 L 

Fig. 1. Eigenvalues tending to be dense in R~ 



5. Uncle Hamiltonians of injective spin chains 

As we have seen in Sec. |3.4[ one can also construct interesting uncle Hamiltonians 
for injective MPS, if one chooses a non-canonical MPS representation of the state. 
To this end, let us start from an MPS \M{A)) with injective tensor A, and let 




Then, \M{C)) = \M{A)) denotes another MPS description of the same state. We 
can now consider a perturbation C* + f ( p p ) and construct the corresponding 
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uncle Hamiltonian H' . Note that not any perturbation would lead to the same 
type of result. The tensor C has additional symmetries since both diagonal 
blocks are the same, and therefore some symmetries are also needed for the 
perturbation. 

The following is the result for the second type of perturbation. 

Theorem 9. Let A he the injective tensor description of a given MPS, and let 
us consider this MPS as described by the non-injective tensor 

A 
A 



Given a perturbation tensor C — ^ p J such that {A R) is injective the ground 
space of the uncle Hamiltonian H' for finite chains is spanned by 



A- A 



and . 

posR 



(32) 



This uncle Hamiltonian is local, frustration free, and gapless. The spectra of 
the finite chains tend to be dense in M"*" . In the thermodynamic limit the ground 
space collapses to a one- dimensional space, which is exactly the ground space 
of the thermodynamic limit of the parent Hamiltonian, and the spectrum is the 

whole positive real line . 

One can check this result following essentially the steps from the preceeding 
sections. In a chain with length A^, one can also consider as low energy states 



the states with momentum \(fik) = X] ■ e' 



2jkm/N 



10), with 



10) 



A 



A- R 



A 



A 



(33) 



The calculations related to states with momentum can be reduced consider- 
ably, since {(pk\Vk) = N{C\^Pk) and {ipk\H'\ipk) = N{C,l\H'\^pk). 
When N is large enough (and odd), we have the approximation 



N 



N 



R n 



R 



A 



E 



A 



R 



A 



R 



. e^'^E^-l 



R 



R 



A 



E 



A 



A 



R 




A 
R 



where 



A 



. The expression multiplying N is either divergent or con- 



vergent to a constant different from for almost every R. Therefore, {ipk\Vk) = 
0{N). 
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The value of {ifk\H'\Lpk) can also be approximated as 

{ifk\H'\^k)^N{(:i\H'\^k) - 




= 0{l/N) 



J 

Hence, there can be found low energy states from this family of states. They 
are orthogonal to the ground space, and can be used to follow the steps in the 
previous sections in order to prove that the uncle Hamiltonian is gapless and 
has spectrum K"'" for most perturbations also for injective MPSs. 



6. Conclusions 



In this work we have shown how to construct new Hamiltonians for Matrix 
Product States, which we called 'uncle' Hamiltonians. These uncle Hamiltonians 
share some of the properties of the parent Hamiltonian (frustration freeness, 
ground space) but have a completely different spectrum. Instead of having an 
energy gap above the zero energy space, they are gapless and their spectrum 
is [0,00). We have shown how these uncle Hamiltonians are obtained by doing 
linear perturbations on the matrices defining the MPS, considering the parent 
Hamiltonian of the perturbed MPS and then making the perturbation tend to 
zero. When the MPS has a gauge invariance, the perturbations leading to an 
uncle Hamiltonian are essentially all that are not compatible with the gauge 
symmetry. This distinction between good directions for the perturbation (those 
preserving the gauge symmetry for which we recover the original gapped parent 
Hamiltonian) and had directions (those incompatible with the symmetry and 
then leading to an uncle) seems to be of upmost importance in the analysis 
of the robustness of topological quantum phases in 2D [5112^. Notice that the 
topological character of a PEPS -the 2D generalization of an MPS- is given 
exactly by the existence of such a gauge symmetry |24j . 
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A. Proof of the intersection and closure properties 



Proof (Intersection property, Lemma^. We start by proving S'2(8iC''nC''(X)S'2 = 
S3. The proof will straightforwardly generalize to Sk <E) (E) Sk = "Sfc+i, 
k>2. 

Let us first show that S'2 C n «) S'2 D S'3 . To this end, let \ip) e S3, i.e., 
there exist X, Y, Z, and W such that 



— ^ 1—^ 1—^ 



r B - B 



B 



pos R 



R B 



E 

pos L 



B 



L - An 



w 



(34) 

where the sums run over the three possible positions of R and L, respectively. If 
we now define 



A 




+ 


R 






1 

X' 


X 






Z 









(35) 



we have that 



A- A 



X 



r A- A 



A 



A 



X' 



e A2 



and similarly 



B - B - B 



Y 



B - B - L 



-W 



e Bo 



E 

pos R^3 



A 



e H2 



and 

pos L^3 



[b}{^ 



■w- 



showing that G ^2 (8 C''. Similarly, one can show that G (E) S2, proving 
that 52 n ® 5*2 D S3. 
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Let US now show that conversely, 5*2 n C"^ ® S'2 C 5*3. To this end, let 
\f) e 5*2 (8) n (g) ^2, i.e., there exist tensors X, Y, etc., such that 



A- A- 



X 



B - B 



A- R 





- A- A- 






B B 


X' 




+ 


Y' 





- L 



L 



Z' 



W 



w 



(36) 



71 1 — 1 1 — I 



w 



B 



We want to show that \ip) is of the form p4| ), i.e., we need to show that X has 
a decomposition such as in (35 1, and so on. To this end, we will make heavy use 
of Lemma [2] In particular, injectivity of {f^ g) implies the existence of a tensor 
left-inverse to R, which at the same time annihilates any of the other tensors 
A, B, and L, as well as the existence of a vector \b) satisfying condition 5 in 

n ® I — 1 _^ r~ 

for rSn = - Id-. We now apply i? to the second site in Eq. (36), 



Lemma 



B 



which cancels all terms except one at each side of (36). Additionally, we project 
the third site onto \b) and obtain 



Z' 



B 



Z' 



(37) 



By defining - Zb 



, we therefore obtain that 



Z' 



(38) 



Similarly, we can apply yl ^ at the second site and a vector \a) verifying condition 

1 @ I — 

5 in Lemma 2 such that rri = - Id 



at the first site to see that 



1 




B 


z 










Z', 



(39) 
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where Z'^ is defined correspondingly. Finally, we can project Eq. (39 onto |5) to 
infer that 

— — (40) 



Corresponding expressions for the f orm of W and W can be derived using L^^ . 

Now let us return to the identity ( 36 1 and apply A^^ to the second site, which 
yields 



A J 



X 



w 



X' 



L R 



Using the analogue of Eqs. (38-40) for W ^ this is equivalent to 



A 



X 



X' 



L z, 



Now, we apply |a) to the first site and obtain 



R 



XI 



where we have defined X '^ ac cordingly. 

Combining Eqs. ([SOl), (iol), and (iSl), we find that 



A 



A 



X 



A 



R 



R 



B 



(41) 



A 



A 



A 



X', 



A 



A 



R 



E 

pos R^3 



A 



R 



B 



which shows that the l.h.s., which is half of the terms in Eq. (34), is contained 



in A3 + -Rj. In the same way it can be shown that the sum of other three terms 
in Eq. (34) is contained in + L3, which proves that Si®'C'^ ® Si C S3. 

The proof that Sk <E) C''' D ^ Sk C Sk+i can be carried out in the same 
fashion, using that the tensors 



^A^ 
E HA 

pos R 



A 



B - B 



R 



B 



B 



E H B 

pos L 



A 



have inverses, since injectivity of ( ^ § ) implies injectivity of 
this can be proven in analogy to Lemma [6j 



(42) 
□ 
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Proof (Closure property, Lemmal^. It is clear that 



Sich n S'right 3 span r A - A 



To show the converse, let |(^) G 5icft n S^-ig 



ight- 



An,r B - B 



B 



1^) = Q 



A 



X' 



B 



Y' 



B 



E 

posR 



R 



B 



Z' 



posL 





1 

- L - 


1 

- A 




*1 




W 









A- X 



Y 



- B 



- B 



pos R, 



A 



R 



B 



pos L 



By applying the inverse tensor corresponding to ^ 

pos R 

2, . . . , A^, we get 





-W 


- B 


- L 


- A- 


- A- 


1 

R 


1 

B 




1 

B -a 



at sites 



Z' 



Now let \a) be such that - A - = - Id- and - B - = (Lemma 2): Projecting 
the first site onto |o) yields -| Z'[- = 0; with a corresponding \b), we find that 

and 



Z - = 0. In the same way, we can prove that -W 

l\T/~nTr tTT/-\ /-trit-t n-nnlTr A 1 -i- /~\ rill on4-/^C' i- 1~\ -R n /-I ^-Vin 





w 




0. 






X 






X' 




Id- 



B ^ to obtain 



Y 



- Y' - = -| Id|- (a similar proof can be found in [21]), 
showing that S'lcfTn 5,.ight C span{|M(A)), \M{B))}. □ 



B. The thermodynamic limit Hamiltonian via the 
GNS-representation 

The way to study the thermodynamic limit of finitely correlated spin chains is 
through a representation of the algebra of local observables [TB] . 

Given a C*-algebra A, and a state uj : A C defined on this algebra, there 
exists a (essentially unique) ^-representation tt of this algebra over a Hilbert 
space H, and a distinguised cyclic vector ^ in this Hilbert space, such that 



u;(A) = (e|7r(A)0, yAeA. 



(43) 



The GNS-representation [B] fulfils this condition. In order to construct it, a 
Hilbert space structure must be introduced via the state lu: {A\B) = uj{A*B) 
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is the (possibly singular) inner product. The quotient of A by the subspace of 
elements such that = is a pre-Hilbert space, which just needs to be 

completed to be the Hilbcrt space H needed for the representation. 

The class to which the identity in A belongs can be taken as the distinguised 
vector ^. 

The algebra we are dealing with is the algebra of local observables over an 
infinite spin chain: 



' I ® I (g) A » Ai+1 



' A- «) I (g) • • • , 



(44) 



i<j 



where each Ak denotes the local algebra of observables at the respective site k. 

Since the dimension at each site is always the same in translationally invariant 

MPSs, this local algebra is the same for every site. 

In the thermodynamic limit setting, it is natural to refer to the limit of the un- 

_L 

B 



normalized states \M{A)) = r A - A 



A 



and \M{B)) = r B 



_L 



B 



as UJA 



A - A- - A- and CJB 



B - B 



B 



., since for any 



local observable O we have Wyi(O) as the limit 



A 



A 



A- A- A 



I I O 



O" I I 



A- A- A 



A- A- A 



because the many copies of the operator 




A 



A- A 
I I O 



A 



A 



A 



O" I I 



A- A- A 



A- A- A 



A 



A 



tend to 




(and a similar 



expression can be used for ujb)- This is exactly the finitely correlated picture of 
[5] which allows to define properly MPS in the thermodynamic limit. 

Let us note that uja and lob are normalized, since the result of the above 
picture using the identity observable -that is, computing the norm- is exactly 
tr{AA) (respectively tr^As)), which was taken to be 1 when we defined the way 
to choose the matrices for an injective MPS. 

For the GNS-representation, we can take uja as the state to which the repre- 
sentation is associatecQ In the first place we have to construct the quotient of 
the algebra A of local observables by the ideal of those observables O such that 
uja{0*0) = 0: 



A- A- A 



A- A- A 



O" 



O' 



O" 



A- A- A 



A - A - A 



^ The clioice of either ujj^ or ujg is irrelevant in this case, since they play a similar role in 
the parent Hamiltonian. However, for general Hamiltonians, the spectrum may depend on the 
ground space taken for the representation. 
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We can see this as the norm of the 'vector' 



O" 



A- A- A 



A 



A- A- 



and the representatives of the equivalence classes in the quotient will have this 
form. Since any operator O can be considered, the Hilbert space of the repre- 
sentation will be seen as the completion of 



S 



where Si 



span ■ 



A 



A 



M 



A 



(45) 

where the corresponding quotient should be taken into account (different tensors 
can lead to the same states). 

The local observables will be represented in B{7i) by themselves tensored 
with the identity, acting on this S = %. 

Since, given any element of 5, only finitely many local Hamiltonians do not 
anihilate on it, the global uncle Hamiltonian is well defined on 5*. If we call H'^ 
the representation of this global Hamiltonian acting on 5, H'^ is densely defined 
on S. 

For H'^ to have a unique self-adjoint extension, it suffices to have a dense set 
of analytic elements in its domain j_^T]. A vector \lp) is analytic if there exists 
some r > such that 

E'^l|ff"(l¥'))ll <oo. 



If \(p) belongs to S it must be in some S-mm- Recall that, for every N, 
\\H\s.^A < 27V + 2 and H{S.n,n) C S-n-i,n+i, and therefore \\H"{\<p))\\ < 
-^fe=i(2M + 2k)\\\ifi)\\. From this we have 



|77"(|<^))||<r"i7, 



2M + 2k 



< r"7T^^i(2M -t- 2) = r"(2Af + 2)' 



which is summable for r < 1/(2M -I- 2). 

Thus every vector in S is an analytic vector for H'^, which is essentially self- 
adjoint, and therefore we can keep calling i?^ its unique self-adjoint extension, 
whose spectrum we must study. 



C. SPECTRAL REPRESENTATION THEOREM FOR 
UNBOUNDED OPERATORS 

The uncle Hamiltonian whose spectrum we want to study is an unbounded self- 
adjoint operator. Therefore its residual spectrum is empty, and the whole spec- 
trum is real. Moreover, since it is a positive operator the spectrum is contained 
in M+. 

The elements in the union of both the continuous and the point spectra 
of an operator A (called respectively (Tc{A) and ap{A)), which in this case is 
the whole spectrum, can be characterized in this way: A S crp{A) U (Jc{A) iff 
there exists a sequence of elements with norm one {\f\,k)}k G ^ such that 
11(^4 — XI)\ipx,k}\\ — ^ 0. Such sequences are called Weyl sequences. In the case 
that A G (7p(A), a constant sequence exists verifying this condition -always the 
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eigenvector. In the case A G (Tc(^) the elements of the Weyl sequences can be 
thought of as 'almost eigenvectors', or 'approximated eigenvectors'. 

Besides the projectors onto the eigenspaces, a spectral measure with some 
associated projectors can be defined. 

Definition 7. JS^ If X is a set, Q is a a-algebra of subsets of X , and H is a 

Hilbert space, a spectral measure for {X, fl, H) is a function E : fl ^ B{'H) such 
that: 

a) for each A in f2, E{A) is a projection; 

b) £(0) = and E{X) = I; 

c) E{Ai n A2) = E{Ai)E{A2) for Ai,A2ef2 (and therefore every pair of 
such projectors commute, since intersection of sets is commutative); 

d) if {An}'^:^i are pair wise disjoint sets from Q , then Ei\J^^^ An) —Yl^=i-^{^ri)- 

And the spectral representation theorem for unbounded operators states: 

Theorem 10. 16^ Spectral theorem (for unbounded operators). If N is a normal 
operator on %, then there is a unique spectral measure E defined on the Borel 
subsets of C such that: 

a) N = J zdE{z) 

b) E{A) ^OifAn a{N) = 

c) if U is an open subset of C and U D <y{N) 7^ 0, then E{U) 7^ 

The spectral representation theorem for bounded operators is completely 
analogous. Using the right one we can proof that the uncle Hamiltonians are 
gapless in the thermodynamic limit - one can sometimes restrict the uncle Hamil- 
tonian to a subspace in which it is bounded. 

With these tools we are able to show that the uncle Hamiltonians are generally 
gapless in the thermodynamic limit. 

Proof: H'^ is gapless. 

The representation of Hp with respect to a ground state w, say iJ^, can be 
uniquely extended to a self-adjoint -and therefore normal- unbounded operator 
[S], with positive spectrum. 

With the spectral theorem above, it can be proven that if we have a unitary 
vector \ip) £ S with ((^jiJ^jcp) = a which is orthogonal to the ground space of 
F^,then iO,a]na{H'J^(l). 

Let us suppose that (t{H'^) C {0} U (a, 00). In such a case the norm of 
would be 

{^W) = {V\ I dE{zM = 1 (46) 

J (a. 00) 

because \(f) is orthogonal to the ground space, and we would have that 
a^{^\H'M^{^\ I zd£;(z)|^) '^^^"=^"^^' 

= (^1 ( / zdE[z) + f zdE{z) \ Iv?) "(^^)^0'"1=» I zdE{z)\ip) > 

J (a, 00) / J (a, 00) 

> {^\ I adE{z)\^) — a{'^\(p) — a. 

J (a, 00) 
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Hence if such |<^) exists there must be some part of the spectrum lying in 



This was the case of the low energy states found for the uncle Hamiltonian 

for the GHZ state, since they were orthogonal to the ground space. However, for 
a general MPS, the low energy states found are not orthogonal to the ground 
space. Therefore we cannot state that E{{0}) in the spectral measure plays no 
role when looking at these states, but we can say that it will be negligible. 

Since the family of states |(/)„) tends to be orthogonal to the ground space, 
we have that (0„|ii^({O})|0„) tends also to 0. We can consider |^„) = — 
E{{0})\(j)n), with norm tending to 1 and orthogonal to lva- 

The only difference is that from the fact that the energies of these states 
arc close to a we cannot infer directly that an element from (0, a] lies in the 
spectrum, but we can prove that some element from (0, r] does for every r > a 
and, therefore, also an element from (0, a] does. If this were not the case, we 
would have that 



(0,a]. 



a 



|j/.n)eker(H-^) 



Jr+ 




which contradicts the hypothesis r > a. 



□ 



